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C^. Elastic charge-exchange in relativistic heavy ion collisions is responsible for 



the non-disruptive change of the charge state of the nuclei. We show that it 
can be reliably calculated within the eikonal approximation for the reaction 

m 

'^ . part. The formalism is applied to the charge-pickup cross sections of 158 

^^ . GeV/nucleon Pb projectiles on several targets. The relative contributions 

o. 

of pion- and rho-exchange are determined, using a single-particle model for 
the internal structure of the nuclei. The calculated cross sections are non- 
negligible for heavy targets. It is shown that these cross sections can be 
^ . useful to obtain information on the Gamow- Teller transition strengths of the 



nuclei. 
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I. INTRODUCTION 

Heavy-ion charge-pickup reactions can help us to access useful information on beta- 
decay transition strengths. At low energy nuclear collisions, it is well known that charge- 
exchange is accomplished by the proton exchange between the nuclei |l[]. At intermediate 
energies ( ~100 MeV/nucleon) nucleon-exchange competes with the charged meson exchange 
0. At relativistic energies ( > 1 GeV/nucleon) , the charge stripping, AZ = —1, cross 
sections will have a substantial contribution from proton removal |^, but the charge pickup, 
AZ = +1, cross section will be solely due to the charged meson exchange, mainly vr^- 
and p^-exchange IQ. It is virtually impossible that a proton-pickup will occur at these 
energies. The strong absorption of heavy ions selects large impact parameters and therefore 
emphasizes the longest range part of the charge exchange force. Since the reaction is very 
peripheral, one expects that the charge-exchange process is practically determined by the 
participation of the valence nucleons. Therefore, these reactions should be a probe of charge 
exchange in a nuclear environment. A series of studies of this type has been performed 
theoretically and experimentally [§]-§]. 

The charge-pickup cross sections for relativistic heavy ion collisions have recently been 
studied at CERN 0-Q. Theoretically, little is known about these reactions. It is the aim 
of this paper to develop a simple description of charge-pickup reactions in relativistic heavy 
ion collisions in terms of what we believe to be the most important ingredients, namely 
the microscopic vr- and p-exchange potentials. An eikonal approach developed in ref. @] for 
charge-exchange in nucleus-nucleus scattering is used. Simple expressions are found which 
can be useful for estimation purposes in the planning of future experiments. 

One basic assumption of our work is that the dominant process is pion and rho exchange 
between a projectile and a target nucleon. However, other mechanisms could also be of equal 
importance, i.e. N + N — > N + A followed by absorption. There has been a quite large 
number of measurements on spin-isospin excitations in charge-exchange ( ^He, t) experi- 
ments at SATURNE that show the physics of these kind of processes are quite complicated. 



We will not study these questions in the present paper, leaving this for a future work. But, 



we refer to some of the SATURNE papers p!3|-p!6|] for the interested reader 



In section 2 we describe the reaction mechanism and the structure ingredients of our 
calculation. In section 3 the dependence of the cross sections on the number of participants 
surface nucleons is studied in details and an application is done for the charge-pickup cross 
sections for 158 GeV/nucleon Pb projectiles on several targets. In section 4 the proportion- 
ality of the measured cross sections to the Gamow- Teller strengths is shown. A comparison 
to the results of section 3 shows that the reaction mechanism does not dependent on the 
structure model. In section 5 we present our conclusions. 

II. RELATIVISTIC HEAVY ION CHARGE-EXCHANGE REACTIONS 

Following ref. [Q, the total cross section for charge-pickup cross sections in high energy 
heavy ion collisions is given by 

a = 271 bV{b)db, (2.1) 

Jo 

where V is the probability of one-boson-exchange at the impact parameter b and is given by 
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V{b) = i^-^^^ ) (2jp + l)-\2jT + 1)-^ exp{-2 Jm x{b)} E 



M(m, u, b) 
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(2.2) 



where Im x{b) is the imaginary part of the eikonal phase. In high energy collisions the 
phase x{b) will be predominantly imaginary and can be constructed from the t-matrix for 
nucleon-nucleon scattering fi% . 

The matrix element M(m, u, b) carries information on the nuclear structure, and is given 
by 

/■oo /'27r 

M(m, z/, fe) = / dqqJ^{qb) / rf0„ e"^'^'^' M{m, q) , (2.3) 

JO JO 

with 

M{m, q) =< $;(rp) ^]{vT)\e-'^-''^ ^(q) e'^'''^ \^^ {r p) ^J{vt) > , (2.4) 



where ^ fh) are the internal wavefunctions of the projectile (P) and target (T) in the initial 
(i) and final (f) states, respectively. 

In eq. ( p.2|) , m = {rrij', mip, mp, m'p) is the set of angular momentum quantum numbers 
of the projectile and target wavefunction. m is measured along the beam axis, and the sub 
indexes T and P refer to the target and projectile, respectively. The probability is obtained 
by an average of initial spins and a sum over final spins {mpijy^ m'p,rp^). The interaction 
potential responsible for the charge-exchange between the nuclei is given by V{q), where q 
is the Fourier transform of the microscopic charge-exchange interaction. For more details, 



see 



The equations (|2.1|) and ( p.2|) are the basic results of the eikonal approach to the descrip- 
tion of heavy-ion charge-exchange reactions at relativistic energies. They can also be used 
for the calculation of the excitation of A particles in nucleus-nucleus peripheral collisions. 
The essential quantity to proceed further is the matrix element given by eq. (|2.4| ) which 
is needed to calculate the impact-parameter-dependent amplitude M(m, u, b) through eq. 
(p.3|). The magnitude of this amplitude decreases with the decreasing overlap between the 
nuclei, i.e. with the impact parameter b. At small impact parameters the strong absorption 
will reduce the charge-exchange probability. Therefore, we expect that the probability given 
by eq. ( p.2| ) is peaked at the grazing impact parameter. 

The pion+rho exchange potential, modified so that the zero range force is corrected for 
an extended source function |]l8l, can be written as 



^(q) = K(q) + Vp(q) 



^(q)(o'i ■ q)(o'2 ■ q) + w(q) (o-i • ^2) 



(ri ■ Ta) 



(2.5) 



where 



^;(q) = vl'"'{q) + <"^(q) , and w{q) = <"*(q) + ^ <"*(q) + <"'(q) + <"'(q) 



(2.6) 



and 



v'n^) 



'j-ji 



q^ 



ml + q' 



2 ' 



<"^(q) 



Jn 



q^ 



mi + q^ 



(2.7) 



wT\^) = -\j. 
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<"*(q) = -Ij, 
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3 " m^ + q^ ' p \ '■' 2 z' ^ _|_ q: 



(2.8) 



<"^(q) = ^ ^^ 3i^ ' <"^(q) = -^ ^P 3^ • (2-9) 



with the parameters g'^ = 1/3, (7|, = 2/3, and E, = 0.4 ]TP| . 

The values of the couphng constants 7,^ and Jp in nuclear units are given by 

J. = 4 = /.^ -^^ ^ 400 Me\/ /m^ 



/p ^^2 (M\ 



J, = ^ = /p 7:^-b - 790 MeV /m^ (2.10) 



where jI/Atx = 0.08 {f'^/4:n = 4.85), m^c^ = 145 MeV, and nipC^ = 770 MeV. 

Turning off the terms w^*^"*, or v^^";'^ and w^^",^, allows us to study the contributions from 
the central and the tensor interaction, and from vr- and p-exchange, respectively. 

Using eq. (|2.5|) , and single-particle wavefunctions, (pjim, it is straightforward to show 
that eq. ( p.4|) reduces to @] (note that here the indexes vr and u denote the proton and the 
neutron wavefunctions, respectively) 

p 
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(2.11) 



Expanding e*'*'" into multipoles we can write 
<05:L'kpe*''-n051>=47rE^'>^;Af(q) <4L'\jiiqr)YjMir)a,\<l,fl> . (2.12) 

IM 

Since ji{qr) Y/Af (i") is an irreducible tensor 

a^ jj{qr) YiM{r) = Yl {nMfi\I'M') ^^a/' , (2.13) 

I'M' 

where "^pm' is also an irreducible tensor. Therefore, 
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a^Jiiqr)YjM{r)\cl>%i>= ^ (J1M^|J'M') < <P^;>^,\<ilpM'\^ye'm> 

I'M' 

J2 {IlMfi\I'M') {jI'mM'lfm') < (f)^f\\^r\\(f)f^ > . (2.14) 



I'M' 



The eqs. ( p .111 - 123^ ) allows one to calculate the charge-exchange between single-particle 



orbitals. The quantity needed is the reduced matrix element < (f)JJ \ \ [ji{qr) a^Yj] 



lM 



>. 



These are calculated in textbooks of nuclear structure (see, e.g., pO[). If several orbitals 
contribute to the process, the respective amplitudes can be added and further on averaged 
in the cross sections. 



III. SINGLE-PARTICLE MATRIX-ELEMENTS AND SURFACE 

APPROXIMATION 

The reduced matrix elements are calculated in the single-particle model for neutrons and 
protons. Using eq. (A. 2. 24) of ref. [^ one finds 
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< 0}' 1 1 [Jiiqr] cy ® Yi]i' 1 1 0f >= - (-1) 



.ii+ii'+j'-l/2 
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^{(2/ +1) //(/+!)} 



1 



\ 



x^+l/2-j 



\{-ir'^'-'m'{22') 



Tijy {q) 
(3.1) 



where 



and 



^/,y (g) = l^ Rf\r)R^\T) ji{qr) r^ dr , 



Vi' Uf) = I {(2j + 1) + (-1)^-*-^'-'' (2j' + 1)} 



(3.2) 



(3.3) 



In eq. ( p.2|) , Rj{r) and Rp{r) are the proton and neutron single particle radial wavefunc- 
tions, respectively. 
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The (3x2) and (2x1) arrays in the equation are matrices and matrix multiphcation 
is imphed. The resulting matrix A^n is a (3x1) array with An, A21 and A^i the values 
of the reduced matrix elements when /' = / + 1, /, and / — 1, respectively. The factor 
1 + (— 1) /2 vanishes unless i + i' + I is even, that is unless parity is conserved. 



The proton and neutron wavefunctions which are needed for the integral eq. ( ^.2] ) must 
carefully account for the Pauli blocking in the final state. Moreover, the pion, or p, is 
readily absorbed in the nuclear surface and only the nucleons in the last shells and close 
to the nuclear surface will contribute to the process. The surface approximation consists in 
using 

47r (Rfir)R^':\r)) = N^ p (r) (3.4) 

where () means the average over the A^^ participant surface nucleons, and is p (r) the 
nuclear density. To simplify, we also use i = £' = 0, so that only the spin of the nucleons 
are considered in the following developments. This will limit J' = 1, and 1 = 0. Thus, the 
relevant matrix element is 

< 0S;) 1 1 lUqr) a ® Fo]i 1 1 0S/^2 >= -^ vf -^0 (?) • (3.5) 

The Woods-Saxon distribution with central density po, radius Rq, and diffusiveness a, 
gives a good description of the densities of the nuclei involved in our calculation. However, 
this distribution is very well described by the convolution of a hard sphere and an Yukawa 
function p3[. In this case, J-'q (g), can be calculated analytically |^ 

1 



•^0 (g) = — — [sin (qRo) - qRo cos (qRo)] 



(3.6) 



Figure 1 compares JFq (g) obtained with the numerical integration with the Woods-Saxon 
density distribution for Al {R^ = 3.07 f m , a = 0.519 fm), Cu {Rq = 4.163 f m , a = 0.606 
fm), Sn (i?o = 5.412 f m , a = 0.560 fm), Au {Rq = 6.43 fm , a = 0. 541 fm) , and Pb 
{Rq = 6.62 f m , a = 0.546 fm) (to simplify the figure, we did not plot the curves for Au and 
Cu). In all cases we used for the Yukawa function parameter in eq. (|3.6|) ay = 0.7 fm. We 



see that the agreement is very good. For carbon, one can use a Gaussian density, with the 
Gaussian parameter a = w| (r^)^ = 2.018 fm^. In this case, 

^0 {q) = n'/'poa'/' exp {-q^a'/A) . (3.7) 

A further simphfication can be obtained for the factor T{b) = exp (— 2/m x{b)) in equa- 



tion ( p. 21) . We use the "tpp" approximation [jl^], which gives 

T{b) =exp|-ajv7V y°° dzj ppir)pT{'R + r)d^r\ (3.8) 

with R = (b, z) . ai\ii^ is the nucleon-nucleon cross section at the corresponding bombarding 
energy, and pp^r) is the projectile (target) matter density distribution. T{b) is known as the 
transparency function. At 158 GeV/nucleon we use a^N = 52 mb. 

The simplest parameterization for the nuclear matter densities are Gaussian functions. 
Assuming 

pp(T) (r) = pPiT) (0) exp [-r"^ / al^T)] (3.9) 

the integrals in eq. ( |3.8| ) can be performed analytically. One gets 

(3.10) 



vr^ aNN Pr(0) Pp(0 ) "Ta? 
(a| + a],) 



T{b) = exp { - -^^^^ r..-^r.^.^^-r-y ^^^ 



(af. + aj,) 



As observed by Karol |22|, for nuclei which densities well described by Woods-Saxon 
distributions, T{b) is very little dependent on the lower values of b and consequently on the 
values of pp^r) (^) for small r's . Only the surface form of the density is relevant. Thus one 
can fit the surface part of the densities by Gaussian functions and use the eq. ( p.lO|) with 



the appropriate fitting parameters pp(0), Pt(0), a^ and ap . If the density distributions is 
given by a Fermi, or Woods-Saxon, function 

p{r) = po{l + exp[(r - R)/{t/AA)]}-' (3.11) 



Karol p2] has shown that T(6) can be well reproduced with Gaussian fits for the nuclear 



densities, if the parameters in the Gaussian distributions (|3.9| ) are given by 
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«' = ^ , p(0) = -poe^/" , (3.12) 

where 

"» = 4^WTFF7T06iP)I ■ * = 4(l„5) = 6.43775... (3.13) 

and A is the mass number. 

In figure 2 we compare the Karol transparency functions (open circles) with the ones 
obtained by a numerical integration of the integral ( |3.8| ) with Woods-Saxon densities (solid 
circles). The agreement is excellent. Thus, the calculation simplifies enormously with the 
use of this approximation. 

In eq. ( p.4|) , only the transverse part of q is needed. Using 

1 1/2 

{i - m)\\ {i + my. 
0, otherwise, (3-14) 



^, r^ ^ , i^(^+m)/2/2£+lV/2[(^-^)!(^ + "^)!] ^ i^S .r ,, 



and 

r27r 



e 







i(— )<A^0 = 27r5^,, (3.15) 



in the matrix elements ( |2.11[ - prT^ ) it is straightforward to show that ( |2.3| ) becomes 



M(m, z/, b) = Co(mp, m'p,mT, ■m'j<)Fo{b) + Cv{mp, m'p,mT, ■m'rp)G„{b), 
where 

Co{mp, m'p,mT, m^ ) = - ^ (-lmp/i|-mp)(-lmT/i|-m^) 

C^{mp, m'p,mT, m'j,) = - ^ (-lmp/i|-mp)(-lmT/u'|-m^)(5^_^/_^ (3.16) 

and 

Fo(6)=/ rfggJo(g&)^(g)^o(?), ^,(6) = / rfg gJ.(g6)i;(g)^o' (?) • (3-17) 

Jo Jo 

In figure 3 we plot the charge-pickup probabilities for Pb (158 GeV/nucleon) + X (tar- 
get), with X = C, Al, Cu, Sn, Au and Pb. The exchange probability is peaked at grazing 
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impact parameters: at low impact parameters the strong absorption makes the probabihty 
small, whereas at large impact parameters it is small because of the short-range of the ex- 
change potentials. The value of the exchange probability at the peak is about 10~^ for Pb 
+ Pb. As in the case studied in ref. Q] the process is dominated by vr-exchange, with only 
a small fraction, of the order of 10%, or less, originating from the p exchange channel. 

The probabilities are divided by the square of the number of participating nucleons. 
As shown, it increases with the target mass number. The cross sections are found out to 
be approximately constant for bombarding energies above 10 GeV/nucleon. The reason is 
simple, since from eq. ( |2.2|) we see that, apart from the factor (1/f ) ~ (1/c) , the only 
energy dependence comes from the total nucleon-nucleon cross section in the absorption 
factor T{h). This is approximately constant at relativistic energies. 

In figure 4 we plot the total charge-pickup cross sections in Pb (158 GeV/nucleon) + X 
(target), with X = C, Al, Cu, Sn, Au and Pb. The dashed line is a guide to the eyes. The 
solid line is explained in next section. The cross section is also divided by the square of the 
number of participating nucleons. It increases steadily with the target mass number. 

In the table below we give the values of cf/N'^ (in mb) for the studied reactions. 



Pb (158 GeV/nuc) + X a/A/'f (in mb) 

C 0.0106 

Al 0.0394 

Cu 0.105 

Sn 0.253 

Au 0.572 

Pb 0.587 



Table I - Total charge-pickup cross sections (per participant nucleon) for Pb (158 
GeV/nucleon) + X (target), with X = C, Al, Cu, Sn, Au and Pb. 
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IV. PROPORTIONALITY TO THE GAMOW-TELLER TRANSITION 

STRENGTHS 



It is clear that the approximation (3^) is very rough and that the assumption of £ = 
nucleons at the surface is not reahstic, specially because the surface nucleons are certainly 
not in an s-wave state. However, these approximations have been useful to extract the main 
features of the reaction mechanism in relativistic vr + p exchange. Some of these features are 
very useful for further theoretical developments. It is clear from figures 2 and 3 that only 
nucleons very close to the surface will contribute to the process. Thus, one can simplify 
these approximations by replacing the nucleon coordinates in eq. (|2.11|) by their surface 
positions, Rp and Rt, respectively. In this case eq. (|2.11| ) reduces to the very compact 
expression 

M = [w{q) + v{q) ] Mgt {P — > P') Mqt {T — > T') e^^'^^ e'"^-^^ (4.1) 



where we also used eq. ( p.l4|) . The Gamow- Teller matrix elements Mgt are expectation 



values of the ar operators. I.e., 

Mgt {A ^ A') = J dh p„, (r) = (A' \ar\ A) . (4.2) 

Inserting eq. ( [4.1| ) into eqs. ( |2.4| ) and (^^.2]), and using the integral 



we get 



/ dcf) exp{~iiy(l)) exp{—iqx cos (f)) =2TxJy{qx) , (4.3) 

V{h) = (^^)' exp{-2/mx(6)} Bgt {P ^ P') Bgt {T ^ T') 

xj:\H{iy,b)\\ (4.4) 



where 



Bgt {A -^ A') = \Mgt {A —. A') \' = \{A' |ar| A)\' (4.5) 
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is the Gamow- Teller transition strength of nucleus A. A sum over final spins and average 
over initial spins is implicit. The function H [u, b) is given by 

/■OO 

H {v, b) = 2n dqq [wiq) + v{q) ] J, {qh) J, [q {Rp + Rt)] . (4.6) 

Jo 

This function is peaked at 6 = Rp + Rt and its width is determined by the range of the 
pion-exchange interaction, i.e., A6 ~ 1 fm. 

The expression (|4.4| ) shows the proportionality between the charge exchange probabili- 
ties and the Gamow- Teller transition strengths in relativistic heavy ion collisions. A similar 



relationship was obtained by Taddeucci et al. ||2J] for (p, n) reactions at degrees for proton 



energies of ~ 100 MeV. For heavy-ion reactions at degrees and bombarding energies of 



~ 100 MeV, a similar relationship was also obtained by Osterfeld et al. ^5|. The validity 
of such a proportionality depends on the factorization of the cross sections into two terms, 
one depending on the nuclear sizes and absorption, and the other on the nuclear structure. 
The solid line in figure 4 displays the cross sections for charge-pickup reactions using the 
proportionality expression ( |4.4| ). The results have been normalized to yield the same mag- 
nitude as (J /Nl for the reaction Pb (158 GeV/nucleon) + Sn. One sees that the agreement 
is excellent, showing that the A-dependence of the process is solely due to geometry factors 
(nuclear transparency and the range of the one-boson-exchange interaction). The oversim- 
plified model of the last section is thus unnecessary, but it was useful to show that one can 



use the proportionality expression ( |4.4] ) to access precious information on the beta-decay 



strengths from charge-pickup reactions with relativistic heavy ions. 

V. CONCLUSIONS 

We described the charge exchange in relativistic heavy ion reactions in terms of vr and 
p exchange and the eikonal approximation. We applied the formalism to the calculation of 
charge-pickup probabilities and cross sections for Pb (158 GeV/nucleon) -|- X (target), with 
X = C, Al, Cu, Sn, Au and Pb. 
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Our model yields probabilities and cross sections which are dependent on the number 
of participating nucleons at the nuclear surface. The cross sections for the process are not 
small, as seen from figure 4. Assuming, e.g., that the number of surface nucleons which can 
contribute to the process is of order of 10 for large systems, one gets cross sections of order 
of 50 mb, or more. 

The calculation is useful to support the proportionality of the measured cross sections 
with the Gamow- Teller matrix transition strengths. These are useful nuclear structure 
information. However, the experiments would have to be able to distinguish these "elastic" 
charge-exchange events from more complicated backgrounds, e.g., charge-exchange with pion 
production. 
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VII. FIGURE CAPTIONS. 



1. The from factor jFg (g) of eq. (|3.6|) , compared to the one obtained by the numerical 
integration using Woods-Saxon density distributions for Al, Cu, Sn, Au, and Pb (to 
simplify the figure, we did not plot the curves for Au and Cu). 

2. Comparison of the Karol transparency functions (open circles) with the ones obtained 



by a numerical integration of the integral |3.8| with Woods-Saxon densities (solid cir- 
cles) . 

3. Charge-pickup probabilities for Pb (158 GeV/nucleon) + X (target), with X = C, Al, 
Cu, Sn, Au and Pb. 

4. Total charge-pickup cross sections for Pb (158 GeV/nucleon) + X (target), with X = 
C, Al, Cu, Sn, Au and Pb. The dashed curve is a guide to the eyes. The solid curve 
represents a calculation of the cross sections divided by the Gamow- Teller transition 
strenghts of the nuclei and normalized to the a/N"^ cross section for Pb -|- Sn. 
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